Quantum state reconstruction with imperfect rotations on an inhomogeneously 

broadened ensemble of qubits 
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We present a method for performing quantum state reconstruction on qubits and qubit registers in 
the presence of decoherence and inhomogeneous broadening. The method assumes only rudimentary 
single qubit rotations as well as knowledge of decoherence and loss mechanisms. We show that full 
state reconstruction is possible even in the case where single qubit rotations may only be performed 
imperfectly. Furthermore we show that for ensemble quantum computing proposals, quantum state 
reconstruction is possible even if the ensemble experiences inhomogeneous broadening and if only 
imperfect qubit manipulations are available during state preparation and reconstruction. 
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I. INTRODUCTION 

The task of reconstructing a physical state from a set of 
measurement data is central to all branches of physics. 
Since full knowledge of the quantum state allows only 
the prediction of a statistical distribution of outcomes 
one must perform numerous measurements on identically 
prepared quantum systems in order to reconstruct the 
state. The abstract notion of the quantum mechanical 
ensemble thus plays a very real part in the business of 
quantum state reconstruction. 

In this paper we shall consider the case of quantum 
state reconstruction within the framework of quantum 
computing. Since the ultimate goal of quantum com- 
puting is to coherently steer the evolution of a quantum 
system towards a specific state wich encodes the answer 
to some computational problem, the final state should be 
the single state that encodes the solution with unit (or 
very nearly unit) probability. Presently quantum state 
reconstruction is an invaluable diagnostic tool for deter- 
mining the fidelity with which basic quantum gates may 
be performed. 

A major challenge in all quantum computing systems is 
that of controlling the interactions between qubits. This 
problem is often tackled by introducing an auxiliary de- 
gree of freedom. Some examples are ion trap quantum 
computing where qubits are coupled via collective vibra- 
tional modes of the trapped ions [l[ and quantum com- 
puting with rare earth ions doped in inorganic crystals 
where qubit-qubit interaction is mediated by a change in 
static dipole moment when a particular ion is promoted 
to an electronically excited state @ ■ Thus although qubit 
dynamics are strictly speaking restricted to the qubit 
Hilbert space, many systems employ a larger space for 
the dynamics with the restriction that at the end of a 
gate operation the effective evolution is described within 
the qubit space. This means that the population of the 
auxiliary level and all coherences between the auxiliary 
level and the qubit levels should be sufficiently small at 
the end of a gate operation and it makes sense to restrict 
oneself to the qubit two-level subsystem when perform- 



ing quantum state reconstruction. 

In quantum computing a lot of effort is put into per- 
forming qubit rotations (or gate operations) with as little 
error as possible. Naturally the same sources of errors 
are present during the reconstructive procedure, however 
it is inappropriate that these errors should detract from 
the gate fidelity. We therefore propose that decoherence 
and loss mechanisms be included in the theoretical ac- 
count for the reconstruction. A further complication oc- 
curs in ensemble quantum computing systems where a 
qubit consists of an entire ensemble of two-level systems 
with almost identical physical characteristics. Qubit ma- 
nipulations in such schemes must be robust to the nat- 
ural variation of system parameters [1, 0, [f|, but will 
inevitably lead to small differences in the evolution felt 
by each member of the ensemble. During readout of 
such systems one obtains in fact the ensemble average 
of the measured quantity and the reconstructive proce- 
dure should take this averaging into account. 

The paper is arranged as follows. In Sec. [IT] we analyze 
the case of single qubit state reconstruction. In Sec. IIIII 
we extend the results of Sec.|IT]to include two-qubit regis- 
ters. In Sec. IIVI we present an application of our method 
to a specific physical quantum computing implementa- 
tion. Sec. [V] discusses the case with correlations between 
the unknown quantum state and reconstructive process 
and Sec. IVII concludes the paper. 



II. SINGLE QUBIT STATE RECONSTRUCTION 

A conceptually simple way of reconstructing the quan- 
tum state of a two-level system is to apply a set of rota- 
tions {Di} and to measure the expectation value of suit- 
able observables in the rotated states. From the mea- 
sured data and knowledge of the rotations performed 
one can infer the full density operator of the unrotated 
state. For instance, in [|| coherent emission of radiation 
from rare-earth ions doped in crystals provides informa- 
tion about quantum coherences and is thus applied to 
reconstruct various quantum states of the ions. For two 
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qubits the same procedure may be employed bearing in 
mind that the two qubits must be rotated differently in 
order to distinguish contributions from the singlet and 
triplet states. The case where both qubits are subject to 
the same rotation is explored by Home et al. in Ref. 
In this case full reconstruction of the quantum state is 
not possible and it is necessary to resort to a numerical 
search. 

We begin by recalling that the density matrix of any 
two-level system may be expanded on the set {a{\ 
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(1) 



with <7]_3 the Pauli matrices and cto the identity. The 
condition Tr(p) = 1 is enforced by setting cq = 1 and 
the further condition that p be positive leaves three real 
parameters fulfilling 
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We can then infer how any initial state transforms under 
the rotation D, . 
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where the bjk follow from Eqs. ([5]). 

We shall find it convenient to represent the state simply 
by 



(7) 



A. Qubit rotations 

In order to proceed with the reconstruction we must 
now perform a set of rotations on identically prepared 
copies of the initial state. When dealing with a sin- 
gle qubit the generic Hamiltonian in the standard basis 
(JO), I 1 )} is gi ven by 



H = |<r • nQ, 



(3) 



where u is the vector of Pauli matrices and h is a unit 
vector. When acting for a time t this Hamitonian effects a 
rotation of the Bloch vector through fit about the vector 
ft. For a perfect rotation of the Bloch vector through 
9i about an equatorial axis ft — (cos(<^j), sm(y>j), 0), the 
rotation operator is given by 



Di = D(<pi,6i) 



cos(6>i/2) ~ie~ lv > sin(0i/2; 
-ie**' sin(0i/2) cos(0i/2) 



.(4) 

Due to the linearity of quantum mechanics it is sufficient 
to know how the Pauli matrices transform under these 
rotations. 

DaiD^ = {cos 2 (0/2) +sin 2 (0/2)cos(2( / 3)} ( Ti 

+ sin 2 (0/2)sin(2^)er 2 

— sin(0) sin(</?)<73 (5a) 

Da 2 D ji = sin 2 (0/2)sin(2(^)cri 

+{cos 2 (0/2) - sin 2 (0/2) cos{2tp)}a 2 

+ sin(0) cos(f)(T 3 (5b) 

Da^D^ = sin(0) sm(tp)ai — sin(0) cos(</?)(72 

COS(0)(T3. (5c) 



From Eq. ^ we may construct the superoperator 
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such that the rotated state may be found as 



< = Bv p , 



(9) 



Let us consider measurements which reveal the popula- 
tion of the |0) and |1) levels of the rotated state, i.e. 



*0 - 2 + C 3 

Pi = 1 - Po = \- 
or equivalently using Eqs. (HI)-© 
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For the special case of Eq. ^ we obtain the explicit 
equations 



- sin a cos ip sin 6 cos ip cos 
sin cos (p — sin 8 cos p — cos 




P 1 
(12) 

It is evident that the two equations are linearly depen- 
dent and thus one must choose at least three sets of rota- 
tion angles (0j, ipi) in order to obtain a coefficient matrix 
of rank 3. The coefficients Cj are then found by sim- 
ply inverting the coefficient matrix of the resulting set of 
equations. 
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B. Relaxation effects 



We now consider the complications of decoherence and 
coupling to auxiliary degrees of freedom. The problem of 
quantum state reconstruction in damped systems is well 
studied, for instance a very general treatment is given in 
Q. Here we shall assume that the decoherence mecha- 
nisms and the coupling to auxiliary degrees of freedom 
are well understood and may thus be modelled accurately 
when describing the full system dynamics which may take 
place in a larger space than the qubit space. The proce- 
dure now follows the same lines as previously. We choose 
a set of rotations {(#, <fi)i} and identify a set of parame- 
ters describing decoherence effects. Taking these effects 
into account we simulate the effect of each rotation on 
each of the basis states 
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by numerical integration of the Lindblad master equation 
p = i[p,H]+£ Iclax {p). (14) 
The function £ rc iax given by 1 

Aclax(p) = -\ J2i C l C ^P + pClC m } + G mPCl, 



(15) 

takes the relevant relaxation mechanisms into account. 
For instance spontaneous transitions from an auxiliary 
state |aux) to the qubit state |0) with decay rate T is 
modelled by a term with 
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Notice in Eq. (| 13|) that the sum now runs from j = 
since some population may be lost to the auxiliary state 
forcing us to relax the normalisation constraint and de- 
scribe the state as a four component vector. From the 
calculated bjk and the measured populations Pi for a sin- 
gle set of rotation angles (6, ip) the equation for the Cj is 
given by 
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Since we are now forced to relax the constraint Tr{p) = 1 
because population may be lost to other levels the two 
equations of Eq. (|17[) are in general linearly independent. 
The minimum set of angles needed to reconstruct such a 
state is two. 



C. Ensemble quantum systems 

As a final level of complication we now consider the 
problem of state reconstruction in the context of ensem- 
ble quantum computing. Explicitly we assume that each 



FIG. 1: The probability that a member of the ensemble 
experiences the error parameter value 8 is given by P(5). 



member of the ensemble reacts slightly differently to the 
external controls and that these differences may be quan- 
tified by some set of parameters <5. For example 8 could 
signify a variation in the level structure brought on by the 
microscopic environment around each ensemble member 
or a variation in the external control field over the spatial 
extent of the ensemble. When performing the qubit rota- 
tions each member of the ensemble experiences a slightly 
different evolution from that of its neighbours so that 
with probability P(S) an ensemble member undergoes 
the evolution 

p -> p R (5) = D(9, if, 5) P D*(0, if, 8). (18) 

The final state is given by the ensemble average 
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p R {8)P{8)d8. 



(19) 



It is a valid point that the effects of 8 would also be felt 
during the preparation of the state p. Means exist, how- 
ever, to eliminate the effects of error Hamiltonians that 
are constant in time and thus to compensate for varia- 
tions in 8 0| . In this section we shall thus assume that 
the ensemble represented by p is uncorrelated with the 
value of 8. In order to proceed with the state reconstruc- 
tion we choose a set of N sample values of the variable 
8 as illustrated in Fig. [1] and assign to each a weight pi 
proportional to the probability P(S) such that 
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(20) 



We then simulate the effects of a rotation D(9, <p, 8) on 
each basis state for each sample point 



N 3 



2 4 / J PibjkQ'k- 

i=l k=0 



(21) 



Again choosing two sets of angles ipi) gives a linear 
system of equations analogous to Eq. (fTT|) for the c$ which 
may be inverted to obtain the initial state v p . 



III. TWO QUBIT STATE RECONSTRUCTION 

In the previous section we asserted that a single qubit 
state can be expanded on the Pauli matrices. The state 
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of a two qubit system may therefore be expanded on the 
tensor product space 



EE' 

i=0 3=0 



(22) 



The matrix representation of Oi (g> aj is constructed as the 
Kronecker product of the two matrices at and (jj. The 
state vector v p now becomes a 16 component vector of 
the Cij. By applying rotations to both qubits we obtain 
the rotated two-qubit state 

PR = {Di®D 2 )p{D\®D\) 

3 3 

= (Di <g> D 2 ) ^2 ^2 CijCTi <g> g j [D\ ® D\) 

3 3 

3 3 3 3 

= EE c y(E 5ifcfjfe ) (E h 3k>°k>) 

i=0 j=Q fc=0 fc'=0 
3 3 3 3 

= ^2 E c v b ^bjk'cr k ® (Tfe'. (23) 

j=0 j=0 fc=0 fc'=0 

As mentioned previously it is important that D\ and D 2 
implement different rotations in order to allow full re- 
construction of the state. As in the single qubit case we 
obtain an expression for the population of the basis state 
\nm) 

3 3 

— EE E E sgn(n,m,k,k')cijbikbjk', 

i=0 ]=Q fc={0,3} k' = {0,3} 

(24) 

where sgn(n, m, k, k') represents the sign of the entry 
(n, m) in the tensor product o% <8 <Jy with k = - 0, 3 and 
k' = 0, 3. Eq. ([M]) is analogous to Eq. flT7|) and the re- 
sults of the previous section now extend naturally to the 
two qubit case. We include the effects of decoherence 
and when working with an ensemble quantum comput- 
ing system we include effects of inhomogeneity over the 
ensemble in 6 choosing again a set of sample points from 
the ensemble probability distribution as shown in Fig. 
[TJ For each set of angles {ip\, 6\, ip 2 , 6 2 ) we obtain four 
equations in the Cy . In order to obtain an invertible sys- 
tem of rank 16 we must thus choose at least four sets 
of angles {^fi,0i,ip 2 ,0 2 ). This approach assumes perfect 
measurement statistics which requires an infinite num- 
ber of measurements in the laboratory and numerical 
precision on the order of the machine precision during 
the integration of the master equation. In practice we 
choose at least five and preferably more sets of angles 
and perform a least squares fit to the obtained data. A 
thorough description of quantum state reconstruction by 
least squares inversion is given in Q . 
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FIG. 2: Left: simplified level diagram of a single rare 
earth ion. Right: qubits are prepared as peaks in a hole 
burnt structure. 



IV. APPLICATION TO QUANTUM 
COMPUTING WITH RARE EARTH IONS 
DOPED IN CRYSTALS 

As an example we shall now apply the methods de- 
veloped in the previous sections to the case of quantum 
state reconstruction in the context of quantum comput- 
ing with rare earth ions doped in inorganic crystals. We 
shall base our study on a rare earth quantum computing 
proposal first developed by Ohlsson et al. Q ■ In this pro- 
posal qubit states are encoded in hyperfine levels of the 
rare earth ions embedded in the crystal. The qubit levels 
are coupled via an electronically excited state accessi- 
ble via optical fields. Due to the large inhomogeneous 
broadening of the excited state qubits may be defined by 
isolated peaks in frequency space as seen in Fig. [21 Each 
qubit thus consists of an ensemble of ions centered at 
a given transition frequency. Qubit-qubit interaction is 
mediated by the dipolc dipole coupling between different 
ions. In rare earth quantum computing arbitrary single 
qubit rotations as described by Eq. (0} may be imple- 
mented using two fields flo and f^i simultaneously [1 01 ] - 
The implementation of single qubit rotations is made ro- 
bust to variations in the excited state shift A, within 
the width of the qubit channel (see Fig. [2]), by sweeping 
the frequency of the laser across the width of the qubit 
in frequency space. Typical parameters for the system 
are fl ~ ~ 1 MHz, A - 100 kHz and (5-10 kHz. 
Recently we have shown that errors due to the inhomo- 
geneous broadening (5 of the hyperfine splitting between 
the qubit levels are in fact dynamically suppressed [ll| 
during single qubit rotations, but they cannot be made 
vanishingly small. 

In the following we simulate reconstruction of the qubit 
state |V>) = l/\/2(|0) + |1)) using three sets of rota- 
tion angles. The inhomogeneous profile is modelled as 
a Lorentzian with variable width and we evaluate the 
fidelity of the reconstruction as 



(25) 



where p is the reconstructed state. In order to evaluate 
the usefulness of our method we first perform the recon- 
struction taking the ensemble averaging into account and 
compare this with the reconstruction assuming 6 = for 
all ensemble members. In the latter case the linear equa- 
tion system of Eq. (|17p is not exactly soluble so a least 
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FIG. 3: (Color online) Fidelity of the quantum state 
reconstruction of a single qubit state as a function of 
the width of the inhomogeneous profile. The dashed 
line is found by averaging correctly over the ensemble 
while the full line is found by ignoring variations over 
the ensemble. 
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FIG. 4: (Color online) Fidelity of the quantum state 
reconstruction of a two qubit state as a function of the 
width of the inhomogeneous profile. The dashed line is 
found by averaging correctly over the ensemble while 
the full line is found by ignoring variations over the 
ensemble. 



squares fit is performed. In Fig. [3] we plot the fidelity 
as a function of the width of the inhomogeneous profile. 
As expected the fidelity drops when the width of the 
Lorentzian increases if we do not properly include the 8 
dependent contributions to the final state. We have also 
simulated the reconstruction of the fully entangled two 
qubit state = 1/2(|00) + |01) + |10) - |11)) using six 
sets of angles. In Fig. 2] we plot the fidelity as a function 
of the width of the inhomogeneous profile. In this case 
the drop in fidelity as the width of the Lorentzian in- 
creases is significant if ensemble averaging is disregarded 
and failure to take this into account during the recon- 
structive procedure may lead to overly pessimistic esti- 
mates of the fidelity of the operation producing the en- 
tangled state. 



V. CORRELATED PERTURBATIONS IN 
PREPARATION AND RECONSTRUCTION 
STAGES 

Consider a single qubit ensemble in the pure spin down 
state = (0, 0,-1). Performing a gate operation cor- 
responds to rotating the Bloch vector as illustrated in 
Fig. [5al Each member of the ensemble follows a distinct 
path along the Bloch sphere and the final state v p is the 
mixed state that results from taking the statistical aver- 
age of all states along the upper blue line. If we now wish 
to reconstruct v p we must perform a number of rotations 
as exemplified in Fig. I5bl but since each ensemble mem- 
ber experiences the same error during the tomographic 
rotation as during the initial rotation, the path followed 
on the Bloch sphere by each member is correlated to the 
path followed during the initial rotation. It is important 
to verify that these correlations do not preclude the iden- 




(a) (b) 

FIG. 5: (Color online) (a): Applying a gate operation 

to a pure state qubit causes ensemble members to 
diverge on the Bloch sphere. Only a few representative 
paths corresponding to the dashed lines in Fig. Q] are 
shown, (b): In the subsequent reconstructive rotations 
the paths followed by different ensemble members are 
correlated with the paths followed during the gate 
rotation. 



tification of the state v p . Intuitively it should be clear 
from Fig.[5]that a similar problem occurs no matter which 
aspect of the Hamiltonian of Eq. ([3]) is affected by 8. Any 
aspect of the external control that varies with 8 will cause 
different ensemble members to follow different paths. 

In the preceding sections we assumed that the state to 
be reconstructed did not posses any correlations between 
the quantum state and 8. Whether dealing with an en- 
semble quantum computing system where a single qubit 
is encoded in an ensemble of physical two-level systems 
or "single instance" systems where a single qubit is en- 
coded in a single physical two-level system, this case is 
well described by the quantum operations formalism [l2| 
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The Ek are linear operators on the qubit Hilbert space 
fulfilling the relation 

J2E k El = l. (27) 

k 

In the ensemble case studied above where ensemble mem- 
bers undergo with probability P(5) the unitary evolution 
U$ we have simply 

E k = ^W)U S . (28) 

This is the case depicted graphically in Fig. [5a| We shall 
now consider the case depicted in Fig. [5b] where the evolu- 
tion of each ensemble member is correlated to the distinct 
evolution experienced by that member during the prepa- 
ration of the state. This case is not contained within the 
quantum operations formalism and so it is not a priori 
clear that state reconstruction is possible in this case. We 
shall now prove that full state reconstruction is in fact 
possible. We can assume that at some time in the past 
the ensemble was in a state possessing no correlations 
with 5. The cumulative effects of all subsequent manipu- 
lations including the gate operation used to prepare the 
state is that with probability P(8) a given ensemble mem- 
ber evolves into v p (S). We proceed by expanding v p (S) 
in powers of the inhomogeneity parameter S 

v p {6) = vf + 6v p V + p v f) + 0(5 3 ). (29) 

We also expand the superoperator B describing the to- 
mographic rotation in powers of 5 

B(S) = B<® + SB™ + 6 2 bW + 0{5 3 ). (30) 

We may think of each ensemble member as being finally 
rotated with probability P(S) into the state 

v*(S) = (B (0) + SB™ + PB&)(vW + SvW + 5 2 v^) 

+S 2 (B^v^ + B^v® + B^v^), (31) 

where we have kept only terms up to 0(5 2 ). When per- 
forming a measurement on the rotated state we therefore 



obtain the ensemble average 

< = J P(6)v*(5)d6 

= B^v^+S(B^p^+B^v^) 

+5 2 ~(BW V W + B&vW + B(%< 3 >). (32) 

If the distribution P(S) is symmetric around the origin 
the first order term averages to zero. However from Eq. 
(|3"2"|) it is evident that v p may be reconstructed to sec- 
ond order simply by expanding the rotation operator to 
the same order. Similarly the state can in fact be found 
to any order as long as the rotation operator is known 
to the same order. In order to determine the individual 
terms in the expansion of the state given in Eq. (|29|) one 
must employ successively more sets of rotation angles. 
If, however we are only interested in the measured en- 
semble average the discussion of Sec. IIIII still holds. In 
this case the expansion of the state is contained implic- 
itly in the averaging over S and four sets of angles is the 
minimum to reconstruct the state. In a single instance 
quantum computing proposal the above arguments show 
that quantum state reconstruction can be performed even 
in systems with intrinsic systematic errors quantified by 
S. 



VI. CONCLUSION 

In conclusion we have described a simple and effec- 
tive method of performing quantum state reconstruction 
in the context of imperfect quantum computing. The 
method can be adapted to include decoherence effects in- 
herent in a given specific quantum computing proposal. 
Moreover we have shown that full quantum state recon- 
struction is possible even in systems that interact with 
environmental degrees of freedom that were also impor- 
tant in the preparation of the state. This is a good exam- 
ple of a process which is not contained within the quan- 
tum operations formalism, and we note that such a sce- 
nario was already mentioned with a rather contrived yet 
illuminating example in 13]. 
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